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ABSTRACT 

We present a possible OVIII X-ray absorption line at z — 0.117 ± 0.001 
which, if confirmed, will be the first one associated with a broad HI Ly/3 (BLB: 
FWHM=160l3Q km s -1 ) absorber. The absorber lies along the line of sight to 
the nearby (z = 0.1372) Seyfert 1 galaxy PKS 0558-504, consistent with being a 
WHIM filament. The X-ray absorber is marginally detected in two independent 
XMM-Newton spectra of PKS 0558-504, a long ~ 600 ks Guest- Observer obser- 
vation and a shorter, ~ 300 ks total, calibration observation, with a combined 
single line statistical significance of 2.8a (2.7a and 1.2<r in the two spectra, re- 
spectively). When fitted with our self-consistent hybrid-photoionization WHIM 
models, the combined XMM-Newton spectrum is consistent with the presence 
of OVIII Ka at z = (0.117 ± 0.001). This model gives best fitting tempera- 
ture and equivalent H column density of the absorber of logT = 6.56lo 17 K, 
and logNf/ = (21.5 ± 0.3)(Z/Zo.oio) _1 cm -2 , and predicts the marginal contribu- 
tion of only two more lines within the XMM-Newton RGS band pass, NelX Ka 
(A = 13.45 A) and FeXVII L (A = 15.02 A), both with equivalent widths well 
within the la sensitivity of the combined XMM-Newton spectrum of PKS 0558- 
504 (EW 1,J < 3 mA). The lack of detection of associated OVI in the archival 
FUSE spectrum of PKS 0558-504, allows us to infer a tighter lower limit on the 
temperature, of logT > 6.52 K (at la). 
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The statistical sigificance of this single X-ray detection is increased by the 
detection of broad and complex HI Ly/3 absorption in archival FUSE spectra 
of PKS 0558-504, at redshifts z = 0.1183 ± 0.0001 consistent with the best- 
fitting redshift of the X-ray absorber. The FUSE spectrum shows a broad 
(FWHM= 160^30 km s _1 ) absorption complex, which we identify as HI Ly/3 
zblb = (0.1183 ± 0.0001). The single line statistical significance of this line is 
4.1a (3.7a if systematics are considered). A possible HI Lya is marginally hinted 
in an archival low-resolution (AA ~ 6 A) IUE spectrum of PKS 0558-504, at a 
redshift of z = (0.119 ±0.001) and with single line significance of 1.7a. Thus, the 
combined significance of the three (XMM- Newton, FUSE, and IUE) independent 
tenative detections, is 5.2a (5.0a if the HI Lya is not considered, and 4.6a if the 
systematics in FUSE are considered). 

The detection of both metal and H lines at a consistent redshift, in this hot 
absorbing system, allows us to speculate on its metallicity. By associating the 
bulk of the X-ray absorber with the BLB line detected in the FUSE spectrum at 
z BLB = 0.1183 ± 0.0001, we obtain a metallicity of 1-4% Solar. 

Although the absorber is only blueshifted by ~ —6000 km s _1 from the sys- 
temic redshift of PKS 0558-504, the identification of the absorbing gas with a 
high velocity nuclear ionized outflow, is unlikely. The physical, chemical and 
dynamical properties of the detected absorber are all quite different from those 
typically found in the Warm Absorber (WA) outflows, commonly detected in 
Seyferts and higher luminosity quasars. WA outflow velocities typically span a 
range of few hundreds to ~ 1 — 2 thousands km s -1 ; WA metallicities, when 
measured, are typically found to be at least Solar; high-ionization WAs are vir- 
tually always found to co-exist with lower-ionization X-ray and UV phases. All 
this strongly suggests that the absorber, if confirmed, is an intervening WHIM 
system. 

Subject headings: WHIM, Absorption Lines 

1. Introduction 

In the present epoch (z< 1 — 2), over half (54 %, Fukugita, 2003) of the baryons are 
missing. They are predicted by Big-Bang Nucleosynthesis (e.g. Kirkman et al., 2003), 
inferred by density fluctuations of the Cosmic Microwawe Background (e.g. Bennet et al., 
2003; Spergel et al., 2007), and seen at z ~ 3 in the "Lya Forest" (e.g. Rauch, 1998; Weinberg 
et al., 1997), but by z < 2 they are unaccounted for in detected stars and gas (Fukugita, 
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2003). According to hydrodynamical simulations for the formation of structures in a A-CDM 
Universe (e.g. Cen & Ostriker, 2006), most, if not all, of these missing baryons should be 
in a barely visible "warm- hot" phase of the filamentary intergalactic medium (the WHIM). 
The WHIM was shock-heated to temperatures of 10 5 — 10 7 K during the continued process of 
collapse and structure formation, and enriched up to Z WHIM = 0.1 — 1 Z by galaxy super- 
winds (GSW, Cen & Ostriker, 2006). These GSWs efficiently removed metals from the cool 
and dense ISM phase in galaxies and spread them into the tenuous (n e ~ 10~ 6 — 10 -4 cm -3 , 
i.e. overdensities 5 ~ 5 — 500 compared to the average density of the Universe) WHIM phase, 
enriching the gas to few percent, and up to tens of percent, of the solar metallicity values 
(e.g. Cen & Ostriker, 2006). WHIM filaments are too tenuous to be detected through their 
bremsstrahlung and line emission with current X-ray instruments (e.g. Yoshikawa et al., 
2003). However, integrated column densities of high- ionization metal ions along a random 
section of one of these filaments could reach values as high as ~ 10 16 cm" 2 (OVII in the soft 
X-Rays) and ~ 10 14 cm -2 (OVI in the Far-Ultraviolet: FUV), imprinting metal absorption 
lines in the FUV and soft X-ray spectra of background sources with equivalent width ranging 
between EW ~ 1-20 mA (the OVII Ka in the X-rays) and EW~ 10-100 mA (the A ~ 1032 
A transition, inthe FUV). 

The most intense of these absorption lines are the OVIi s 2 2s _ 5 .i s 2 2p doublet in the UV 
(A = 1031.926, 1037.617 A), and the CV Ka(r) (A = 40.268 A), OVII Ka(r) (A = 21.602 A), 
CVI Lya (A = 33.74 A) and OVIII Lya (A = 18.97 A), in the X-ray band. Which of these 
lines dominates, depends on the ionization state of the gas, that is mainly on its temperature 
(and at a second order on the gas volume density). Dave and collaborators (2001) showed 
that the baryon temperature distribution in the intergalactic space, peaks at logT ~ 6.6 K, 
and is strongly skewed toward low temperatures: 50% of the baryon in the WHIM are found 
between logT = 5.6 K and logT = 6.7 K at a weighted average temperature of logT = 5.9 
K, while the low (logT = (5 — 5.6) K) and high (logT= (6.6 — 7) K) temperature tails of 
the WHIM temperature distribution, contain 27% and 23% of the gaseous baryons in the 
Universe, respectively, around weighted average teperatures of logT = 5.4 K and logT = 6.7 
(Conciatore et al., in preparation). So, while the vast majority (73%) of the baryons in 
the WHIM should absorb and emit in the X-rays at the wavelengths of the He-like (50%) 
and/or H-like (23%) ions of C and O, a substantial fraction of the WHIM (27%) should be 
detectable in the FUV, through Li-like C and O transitions. 

These theoretical predictions have been confirmed by FUV observations (e.g. Danforth 
k Shull, 2005, 2008; Tripp et al., 2008; Richther, 2006). Danforth & Shull (2005) analyzed 
FUSE data of 31 AGNs with z < 0.15 to search for OVI absorption counterparts to 129 
known Lya absorbers. They found 40 such systems, deriving a dNovi I 'dz{> EWrhresh) 
that agrees strikingly with the latest predictions by Cen & Fang (2006; see their Fig. 2). 



-4- 



However, as pointed out recently by Tripp et al. (2008), the vast majority of the associated 
HI Lja absorption lines are too narrow to be produced in gas with typical WHIM temper- 
atures. Danforth & Shull (2008) proposed that low-ionization metals and narrow HI Lja 
trace mildly photoionized medium (the low- z Lja Forest: ~ 30% of the baryons), while 
narrow-HI/OVI/NV associations trace a multiphase intergalactic medium, with a photoion- 
ized portion of the intervening filament imprinting the narrow-HI absorption and a shock- 
heated part (WHIM) imprinting the OVI/NV on the UV spectra (~ 10% of the baryons). 
Danforth & Shull (2008) estimate ft^j IM = 0.34 % (down to logN y/ > 13.4 cm" 2 ). This 
is ~ 15 % of the "missing mass" (strictly speaking, only a lower limit, given the large un- 
certainties in the ionization correction for the OVI-bearing gas), in good agreement with 
the WHIM baryon fraction predicted to reside in the low-temperature tail of the WHIM 
mass-temperature distribution. 

Detecting the bulk of the 'Missing Baryons' in the "X-Ray Forest" , instead, has proven 
to be extremely difficult. This is because of the unfortunate combination of (a) the still 
limited resolution (R ~ 400 at A = 21.6 A) and the low throughput (Ae// ~ 20 — 40 cm 2 ) 
of the current high-resolution X-ray spectrometers [the XMM-Newton Reflection Grating 
Spectrometer (RGS, den Herder, et al. 2001) and the Chandra Low Energy Transmission 
Grating (LETG: Brinkman et al., 2000)]; (b) the lack of bright {f . 5 -2keV > W 11 erg s^ 1 
cm -2 ) extragalactic point-like targets (e.g. Conciatore et al., in preparation); and, (c) the 
dramatic steepening (Aa>1.5) of the predicted number density of metal WHIM filaments 
per unit redshift at ion column densities >10 15 cm -2 . 

Current evidence is still limited, and highly controversial: in 2005, Nicastro et al. 
(2005a,b) reported the first two > 3a detections of absorption systems identifiable with 
WHIM filaments, at z — 0.011 and z = 0.027, towards the blazar Mkn 421, which was in 
outburst. One of these two systems is at a redshift consistent with that of a known interven- 
ing HI Lja absorber (Shull, Stocke & Penton, 1996) that line, however, is too narrow, and so 
the absorbing gas too cold, to be physically associated with the X-ray metal filament. This 
claim has been questioned in two subsequent papers by Kaastra et al. (2006) and Rasmussen 
et al. (2007), based on: (a) apparent wavelength inconsistencies for some of the lines, which 
called into question their association with a given WHIM system; (b) statistical arguments 
for a significantly lower detection significance for the WHIM identifications than originally 
reported by Nicastro et al. (2005a,b); and (c) the lack of detection of these two systems in 
XMM-Newton RGS calibration spectra of Mkn 421 (but see also Nicastro, Mathur & Elvis, 
2008, for rebuttal arguments). Fang, Canizares & Yao (2007) proposed the detection of an 
OVIII WHIM filament along the line of sight to the blazar PKS 2155-304, once again near 
the redshift of a known narrow HI Lya absorber not due to the same gas as the putative 
X-ray filament, because of the line width. The above reported pieces of evidence for the exis- 
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tence of an "X-ray Forest" at low redshift, have been gathered by exploiting the technique of 
observing random lines of sight toward the brightest possible X-ray sources, with exposures 
long enough to reach extremely high S/N spectra. This technique has the advantage of not 
being biased toward the strongest WHIM systems, and so allows probes of the bulk of the 
WHIM mass distribution. However, the WHIM is predicted, and in the FUV has most likely 
been found (e.g. Stocke et al., 2006b), to correlate strongly with galaxy concentrations in 
the Universe: the densest of such concentrations are also supposed to host the nodes, and 
so the densest and hottest parts, of the WHIM network. 

An alternative technique, therefore, is to target the WHIM search in the X-rays along 
those lines of sight toward bright background X-ray sources that intercept large scale struc- 
ture concentrations. By exploiting this technique Buote et al. (2009), reported a ~ 3a 
detection of a high column density (Novn > 10 16 cm -2 ) OVII Ka absorption line, along the 
line of sight to the blazar H 2356-309, at the redshift of a very large concentration of galax- 
ies in the Sculptor Wall. This technique has the disadvantage of probing only the extreme 
high-temperature/density and low- mass fraction ends (~ 23%) of the WHIM distribution. 

All the above proposed WHIM detections so far, suffer the serious handicap of lacking 
the detection of a clearly associated HI counterpart. The detection of broad HI Lyman lines 
(mostly Lya at A = 1215.67 A - BLA - and Ly/3 at A = 1025.72 A - BLB) is vital for a 
proper assessment of the WHIM metallicity and mass, and require the use of both X-ray and 
FUV facilities. Without HI, metallicity cannot be determined and so a cosmological mass 
density 1 cannot be derived. 

The expected thermal FWHM of an absortion line imprinted from a gas of particles of mass 
m at the equilibrium temperature T, is given by 2\/2Zn2 times the variance of the velocity 
distribution (the Doppler parameter 6), that is FWHM= 2\j2ln2^JkT /m. For HI, this 
implies a broadening of FWHM~ 380 km s _1 at logT = 6.5 . 

The HI fraction, relative to the total H, in gas with temperatures in the logT = (5.8 — 6.3) 
range (where 50% of the WHIM should be, e.g. Dave et al., 2001), spans the interval 
10 -6 ' 2 — 10 -72 . Thus only very shallow (EW(Lya ~ 8 — 80 mA) and broad HI absorption 
lines are expected to be imprinted onto FUV spectra by WHIM filaments, requiring spectra 
with S/N~ 10 — 100 per resolution element to make them detectable. 

Here we present the first tentative detection of a hot (X-ray) WHIM filament along the 
line of sight to the z = 0.1372 Seyfert galaxy PKS 0558-504, that is physically associated with 



For heavier elements, the width scales with the inverse of the root square of the atomic weight. So, for 
example, for O, Ne and Fe, at logT = 6.5, FWHM=95, 85 and 51 km s" 1 . 
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a BLA HI absorber. A full band analysis of the XMM-Newton RGS spectrum of PKS 0558- 
504, is deferred to a companion paper (Papadakis et a!., 2009, submitted). In §2 we present 
the XMM-Newton, FUSE and IUE data that we use in our analysis. §3 and 4 are dedicated 
to the X-ray and FUV data reduction and analysis. In §5 we critically discuss our findings, 
and in §6 we summarize our conclusions. 



2. XMM.- Newton, FUSE and IUE Data of PKS 0558-504 

Between February 2000 and October 2001, PKS 0558-504 was the subject of an extensive 
calibration campaign with XMM-Newton, with a total expsoure of 312 ksec (hereinafter 'CAL 
spectrum'). More recently, in September 2008, PKS 0558-504 was re-observed by the XMM- 
Newton for 618 ksec, as part of an approved cycle 8 Guest Observer program (PI. S. Papdakis; 
hereinafter 'GO spectrum'). In this work we make use of the RGS spectra extracted from 
both data sets. We focus mostly on the analysis of the longer and much higher S/N 2008 (16 
vs 9, at 25 A) GO spectrum. We use the lower S/N 2000-2001 CAL spectrum for consistency 
checks. 

PKS 0558-504 was observed by the FUSE satellite twice: on 1999, December 10 and 
2001, November 7, with net exposure times of 47.2 ksec and 48.3 ksec, respectively. In our 
analysis we co-add the two spectra from the Lif2A channel (with the highest efficiency in the 
1100-1200 A interval) to maximize the S/N and perform our search for intervening BLAs 
and BLBs in this co-added FUSE spectrum @ . 

Finally, PKS 0558-504 has been observed three times with the 'Short- Wavelength' spec- 
trometer of the IUE satellites (the only one covering the wavelength range of interest), on 
1987, September 22, and 1989, November 14-15. In our analysis we use the highest S/N of 
the three spectra (with a net exposure time of 16.2 ksec), the one taken on 1989, November 
15, and check the other two spectra for consistency. 

Table 1 lists the journal of the XMM-Newton, FUSE and IUE observations used in this 
work. For all our spectral fitting we use the fitting package Sherpa, of the Chandra Interactive 
Analysis of Observation (Ciao) software (v. 4.1.2; Freeman, Doe & Siemiginowska, 2001). 



2 We also reduced and analyzed the sum of the two spectra from the LiflB channel, which is sensitive 
(though with about half of the effective area of the Lif2A channel) in a wavelength interval similar to that 
of the Lif2A. However, the LiflB channel is unreliable at > 1140 A (The FUSE Data Handbook, Chapter 7, 
§7.3.2, http://archive.stsci.edu/fuse/DH_Final/Factors_Affecting_FUSE_Data.html), which is the region of 
interest for our analysis. Therefore we decided to consider only the Lif2A channel in our analysis (but see 
also footnote 9, §4). 
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3. The RGS Spectrum of PKS 0558-504 

3.1. RGS Data Reduction and Analysis 

The 2008 GO RGS data of PKS 0558-504 were reduced as in Papadakis et al. (2009), 
and the same standard reduction procedure was applied to the 2000-2001 calibration obser- 
vations of this target. The data were cleaned for periods of high background activity, by 
excluding events taken within time-intervals where the background deviated positively by its 
average level by more than 2a. Finally RGS1 and RGS2 spectra and responses of the single 
observations of the two datasets (CAL and GO) were co-added and averaged by using the 
rgscombine tool. This produced total net exposures of 480 ksec (RGS1) and 466 ksec (RGS2), 
for the GO dataset, and 309 ksec (RGS1) and 283 ksec (RGS2) for the CAL dataset. Despite 
a 30% higher average flux level of PKS 0558-504 during the CAL observations, compared 
to the GO observations, and a relatively small difference of only 65% in the net exposure 
times of the two datasets, the S/N of the GO spectra is significantly (factor of ~ 2) higher 
than that of the CAL spectra. In half RGS resolution element (30 mA), at 0.5 keV the 
co-added RGS1 CAL spectrum has S/N=9, compared with S/N=16 for the GO spectrum. 
This difference in S/N is due to the much higher (more than an order of magnitude) level 
of the average background during all individual CAL observations, compared with the GO 
observations. This is clearly seen in Figure 1, where we plot the combined GO RGS1 20-24 
A background-subtracted source (green filled circles and errorbars) and background (blue 
filled circles and errorbars) spectra of PKS 0558-504, together with the combined CAL RGS1 
20-24 A background-subtracted source (black empty squares and errorbars) and background 
(red empty squares and errorbars) spectra of PKS 0558-504. We therefore decided not to 
further combine the RGS spectra of the two datasets, and to analyze them separately and 
simultaneously. 

In our spectral analysis we bin the RGS spectra two 45 mA bin (about half the nominal 
FWHM resolution of the RGSs). 

3.2. X-Ray Spectral Fitting 

In our companion paper (Papadakis et al., 2009), we present the broad-band 0.4-2 keV 
spectral analysis of the 2008 RGS spectra of PKS 0558-504, and show that the soft X-ray 
continuum of PKS 0558-504 is well modeled by a broken power law with a break energy 
of Eb r k ~ 1.5 keV, attenuated by Galactic ISM absorption. More importantly we exclude 
the presence of a Warm Absorber with typical column densities and ionization parameters 
(Papdakis et al., 2009). This is consistent with the absence of associated, intrinsic, Narrow 
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Absorption Lines (NALs) in the Far Ultraviolet spectrum of PKS 0558-504 (Dunn et al., 
2007). However, Papadakis et al. (2009) also note that weak absorption line-like features 
are detected in the spectrum of PKS 0558-504, in the narrow 20-24 A spectral interval. These 
three lines, and in particular the line at A = 21.18 A are the subject of this study. 

We first explored the 20-24 A region of the GO and CAL RGS1 spectra of PKS 0558-504 
(the RGS2 is blind in this wavelength range, due to the lack of a read-out CCD chip). We 
fitted the data simultaneously, with two local continuum models, each including a power law 
attenuated by neutral absorption with the metallicity set to Anders & Grevesse (1989). While 
such a model gives a statistically acceptable fit for the GO spectrum of PKS 0558-504, it 
leaves broad and systematic residuals in the CAL spectrum. This is because the rgscombine- 
averaged CAL response matrix fails to properly reproduce the depth of the several effective 
area features present in this portion of the RGS1 spectrum. We therefore added a number of 
FWHM~ 0.5 — 1 A broad emission and absorption gaussians to the best fitting continuum 
model of the CAL spectrum of PKS 0558-504, until we reached a statistically acceptable 
fit. Three narrow regions with a deficit of counts are left in the GO and CAL residuals to 
the best fitting continuum models, at A ~ 23.5 A, A ~ 21.6 A and A ~ 21.2 A. To model 
these residuals we added three Gaussians to our best-fitting continuum models and refitted 
the data leaving all Gaussians parameters free to vary independently in the fit to the GO 
spectrum, but linking the line energies and Gaussian widths in fitting the CAL spectrum to 
those of the corresponding parameters of the fitting model of the GO spectrum. The first 
three rows of Table 2 show the best fitting parameters of these three absorption lines. Errors 
on wavelengths and redshifts are assumed to be equal to 1 bin-size, 30 mA, if unresolved by 
the fitting routine). Figure 2 shows the 20-24 A RGS1 CAL and GO data of PKS 0558-504, 
with their best-fitting models superimposed. The first two of these lines, at A = 23.52 ±0.03 
A and A = 21.61 ±0.03 A, are easily identified with the Ka inner shell transition from atomic 
OI in the Interstellar Medium of the Galaxy, and the Ka He-like transition from OVII either 
in our Galaxy halo or in the Local group (e.g. Bregman et al., 2007 and references therein). 
The third line, at A = 21.17 ± 0.04 A, appears marginally resolved (but still consistent with 
an unresolved Gaussian at a la confidence level □ ) in the RGS1 (FWHM<1200 km s" 1 : 
Table 2) and has no obvious identification with transitions at z ~ 0. If red-shifted to the 
systemic redshift of PKS 0558-504 [z = 0.1372, REFs), this line would have a rest-frame 
wavelength of A = 18.62 A, close to the rest-frame wavelength of the OVII K/3 transition. 



3 the projection routine in Sherpa, which we use here to compute uncertainties and upper limits, finds only 
an upper limit for the line FWHM, at a la confidence level. This is because of the high non-gaussianeity 
of the RGS Line Spread Function (e.g. Williams et al., 2005), which makes the effective LSF FWHM 
significantly broader than the nominal Gaussian equivalent LSF FWHM of 60 mA(corresponding to ~ 900 
km s" 1 at 21.2 A). 
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However, in this much stronger (factor ~ 6.4) OVII Ka line □ with EW~ 60 mA 

should be visible at A = 24.55 A. Such a line is not seen in the data (see also Papadakis et 
al., 2009) at a 3a limit of EW< 18 mA[f| . Moreover, the absence of intrinsic OVI absoprtion 
in the FUSE data of PKS 0558-504 (Dunn et al., 2007), makes this interpretation less likely. 

The next possible identification is with an intervening absorber, between the redshift of 
PKS 0558-504 and us. This imposes the condition that the line has a rest-frame wavelength 
in the 18.6-21.2 A spectral interval. The strongest resonant absorption line expected in 
this spectral range is the OVIII Lya transition, with an unresolved doublet at an oscillator 
strength weighted average wavelength < A >= 18.97 A. This would give a redshift of z = 
0.116 ± 0.002 for the intervening absorber. 



3.3. Consistency Test with a Self-Consistent WHIM Spectral Model 

To check whether the presence of an OVUIr-q. WHIM filament at z — 0.116 ± 0.002 
along the line of sight to PKS 0558-504, is consistent with the full band XMM- Newton RGS 
spectra of this target, we used our hybrid ionization (collisional + photoionization by the 
meta-galactic UV and X-ray background at a given redshift) WHIM spectral models (an evo- 
lution of our PHotoinized Absorber Spectral Engine - "PHASE" -, Krongold et al., 2003), to 
simultaneously fit the RGS1 and RGS2 Cal and GO data of PKS 0558-504. The fitting model 
includes our best-fitting continuum model (a combination of power-laws and broad gaussians 
to cure residual calibration uncertainties), two negative gaussians to model the Olxa and 
OVII^q absorptions at z ~ 0, and a WHIM model. In the X-ray band alone one cannot 
constrain the absolute metallicity of the absorber (for which detection of HI is needed), since 
the ionization balance of an hybridly ionized cloud of gas (i.e. the fractional ion abundances), 
parameterized by the temperature of the absorber, at equilibrium, is virtually independent 
on metallicity (see also §5.2). Thus, for fitting purposes only, we froze the WHIM metallicity 
to Solar. The turbulent velocity of the absorber was frozen to 400 km s -1 (about 1/2 of 
the RGS FWHM resolution at 22 A). In a hybrid-ionization WHIM model, photoionization 
by the meta-galactic UV and X-ray photon field, is only a secondary ionization mechanism, 
and its realtive importance, compared to the collisional shock mechanism, depends only on 



4 together with several other strong lines from He- and H-like ions of C, N and Nc 

5 we also caution, however, that although the RGS1 spectra of PKS 0558-504 are free of instrumental 
feature at, exactly, 24.55 A, this spectral region of the RGS1 contains several instrumental features, which 
might reduce, due to calibration uncertainties, the true sensitivity of the RGS1 spectrum at 24.55 A. See 
also §5.2 
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the volume density of the WHIM filament: the lower the density the more important the 
photoionization contribution. In such models the absorber density is thus highly degener- 
ate with the absorber temperature and equivalent H column density and so is difficult to 
constrain. We therefore froze the density of the absorbing gas to a typical expected WHIM 
value of n& = cm~ 3 (overdensity 5 ~ 50 compared to the average baryon density in the 
Universe), corresponding to an ionization parameter (the ratio between the photon density 
at the surface of the absorbing cloud and the baryon density in the cloud) of logU = — 1.9 
at z — 0.1. At such densities (and higher), the photoionization contribution is essentially 
negligible, and the model reduces to a pure collisional ionization model. The remaining 
WHIM parameters, namely the gas temperature and equivalent hydrogen column density, 
and the redshift of the WHIM filament, were left free to vary in the fit. The best-fitting 
WHIM model (Fig. 3, black curve) has logT = 6.56^, logN^ = 19.5 ± 0.3 (for Z = Z Q ) 
and z = 0.117 ±0.001. This model fits the most prominent OVIII Ka line at z — 0.117, and 
predicts only two additional weak absorption lines (EWtl mA): the NeX Ka, at A = 15.02 
A (WHIM filament frame) and a strong (oscillator strength 2.95) FeXVII L line at A = 16.77 
A(Fig. 3, black curve). 

Figure 4, 5 and 6 shows the two portions of the RGS CAL and GO spectra of PKS 0558- 
504, in which the lines predicted by the best-fitting WHIM model lie. The superimposed 
red curve is the best fitting WHIM model convolved with the instrumental responses. The 
OVIII Ka line is detected both in the CAL and GO spectra (Fig. 4). The position of the 
predicted NelX Ka and FeXVII L lines is marked in Fig. 5 (RGS1 and RGS2 GO spectra of 
PKS 0558-504) and Fig. 6 (RGS1 ands RGS2 CAL spectra of PKS 0558-504): the strength 
of these lines is < 1/3 the minimum detectable la EW in the combined CAL+GO spectrum 
(Table 2). We note the hint of the presence of the FeXVII L line at 16.77 A in both the RGS1 
GO (Fig. 5, top panel) and RGS2 CAL (Fig. 6, bottom panel) spectra of PKS 0558-504. 

In Figure 3 we plot the best fitting WHIM model (central, black curve) together with its 
negative (lower, red curve) and positive (upper, blue curve) la temperatures. Both lower and 
upper limits on the temperature are weak. The strongest X-ray lower limit diagnostics in the 
model is the increasing strength of OVII K series absorption lines. However, the sensitivity 
of the RGS1 GO and CAL spectra of PKS 0558-504 at the wavelength of the strongest line 
of this series (the OVII Ka line) is dramatically reduced by the presence of an instrumental 
feature due to a bad pixel in the read out dispersion detectors. On the other end, at higher 
temperature, the limit is set by the decrease of opacity at the wavelengths of the OVIII 
Ka transitions, which is only detected at a combined single-line significance level of 2.8a in 
the data. This can be also seen in Figure 7, where we show the 68%, 90% and 95% logT 
- logN# contours of the best-fitting WHIM model. The parameters are clearly correlated 
and the temperature poorly constrained both at low and high temperatures. However, low 
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temperatures and high column densities are excluded at high confidence levels. 

4. The FUSE and IUE Spectra of PKS 0558-504 

We searched the The Multimission Archive at STscI (MAST) archive | for FAR-UV 
data of PKS 0558-504 and found several data sets, taken with FUSE, HST and IUE satellites 
(Table 1). Our main goal was to search for HI (Lya and Ly/3) and OVI absorption either 
at or near the systemic redshift of PKS 0558-504, or at a redshift consistent with that of 
the putative X-ray WHIM filament, tentatively detected in the RGS spectrum of PKS 0558- 
504. The available archival HST spectrum of PKS 0558-504 was taken with the STIS-G230 
gratings, and so covers a wavelength interval (A = 2760 — 2910 A) not relevant to this 
analysis. The FUSE and IUE spectra, instead, cover two important regions: (a) the HI Ly/3 
and OVI regions (FUSE) and (b) the HI Lya region (IUE). We retrieved these data sets 
from the MAST archive, and analyzed them to search for a secure identification of the X-ray 
absorber. 

FUSE observed PKS 0558-504 twice, in two different programs (P101, Sembach 1999, 
and C149, Prochaska 2002). All FUSE observations after the end of operations were re- 
reduced using the final version of the CalFUSE pipeline (v3.2.3, Dixon et aL, 2007), and we 
used the resultant science calibrated files. We used specific IRAF scripts LI to co-add the 
data from the Lif 2A channel, which covers the wavelength range 1087-1181 A, with a total 
exposure time of 90.3 ks. The spectrum was then binned to two different resolutions: the 
natural resolution of FUSE (3 bins; 11 km/s) and a poorer resolution (20 bins; 70 km/sec) 
with higher S/N per bin, specifically to search for broad, shallow absorption features. 

We use the fitting package Sherpa in Ciao (v. 4.1.2; Freeman, Doe & Siemiginowska, 
2001), to search the FUSE spectrum for: (a) AGN intrinsic absorption between -1000 and 
+1000 km s~ x of the systemic redshift of PKS 0558-504; (b) intervening absorption at a 
redshift consistent with that of the putative X-ray WHIM filament. Specifically, for each 
of the two possible systems, we searched for absorption lines from: (1) the NIII triplet at 
A = 989.80 A, A = 991.51 A, and A = 991.57 A; (2) the OVI doublet at A = 1031.93 A 
and A = 1037.62 A; and (3) the HI Ly/3, at A = 1025.72 A. We confirm the non-detection 
of intrinsic (i.e. nuclear) ionized or neutral absorption, down to EWt30 mA (3a) of Dunn 
et al. (2007). However, we do detect several absorption lines in the wavelength range 1137- 



6 http://archive. stsci.edu/index. html 
7 http://fuse.pha.jhu.edu/analysis/IRAF^cripts.html 
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1155 A (Figure 8). The strongest of these lines are easily identified as Galactic Fell, PII and 
possibly H2 transitions. These lines are all relatively narrow , with typical FWHM~ 50 — 70 
km s _1 , and exhibit two velocity components at v i ~ and v 2 — +200 km s -1 (for the 
low-ionization atomic transitions) □ , and v i ~ and v 2 — +70 km s -1 (for the molecular 
transition: Wakker, 2006): Figure 8. 

However, we also detect a broader absorption complex at the centroid wavelength of 
A = 1147.1 A, and at single-line statistical significances of 4. lex (Figure 8) . We model 
this absorber with a negative Gaussian, and list its best-fitting parameters and equivalent 
width in Table 1 . 

This absorption line is broad (FWHM= 160l 3 o km s _1 ) and cannot be identified with 
any strong Galactic transition. We identify this absorber as a broad HI Ly/3 (BLB) at 
zblb — (0.1183 ± 0.0001). The X-ray redshift of the OVIII Ka absorber is consistent with 
zblb- At the best fitting X-ray temperature of logT=6.56 K, the thermal speed of protons 
is \/2b = 250 km s -1 and so the expected broadening of HI absorption lines is FWHM~ 400 
km s _1 . The width of the zblb HI absorber is marginally consistent, at a 2.5cx level, with 
the expected broadening (Fig. 9). Here we therefore assume that the HI BLB is associated 
with the bulk of the X-ray OVIII absorber. 

There is no detectable OVI or NIB absorption associated with the zblb HI absorber. 
The non-detection of OVI puts a much more stringent lower limit on the temperature of the 
detected WHIM filament. At the best fitting X-ray temperature the relative fraction of OVI, 
£ OVI = 2.0j|- x 10" 4 , predicts (for unsaturated lines) EW(OVI)= 7.6+^° 8 6 mA. From the 



8 With the exception of the high velocity component of the FeII(A1144.94) absorber, at A = 1145.7 A, for 
which we measure FWHM= (100 ± 40). This absorber could be at least partly contaminated by intervening 
HI Lybeta absorption at z = (0.1170 ± 0.0001), i.e. still consistent with the X-ray redshift of the putative 
OVIII intervening absorber, and so possibly representing a second HI component associated with this hot 
intervening WHIM filament. 

9 There are also two known velocity components in high-ionization OVI Galactic (or Galactic halo) ab- 
sorption, along the line of sight to PKS 0558-504, in the velocity ranges —115 — 1-135 km s _1 (Savage, 2003) 
and +210 - +315 km s" 1 (Sembach et al., 2003). 

10 The associated la error here, and anywhere else in the paper, is only statistical. By taking into account 
the vagaries of continuum placement and fixed-pattern noise in the Lif2A FUSE spectrum of PKS 0558-504, 
we estimate an additional 14% systematic la uncertainty. This would reduce the significance of the broad 
line at A = 1147.1 A, from 4.1a to 3.7cr. 

n We note that we this line is also clearly visible in the co-added LiflB spectrum, with a slightly larger, 
but consistent, EW. However, due to the unreliability of the LiflB channel at A^ 1140 A (see footnote 1, 
§2), we conservatively decided not to consider this as a supporting evidence. 
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FUSE data we have EW(0VI) 2 =o.ii8~12 niA at a la confidence level (Table 2). By linearly 
interpolating between expected OVI fractions and EWs, we derive the la upper limit on 
the fraction of OVI corresponding to the measured EW(OVI) upper limit, and from this 
the la lower limit on the temperature of the absorber. This gives: logT> 6.52. Thus by 
combining the UV and X-ray constraints on the temperature of the WHIM filament, one 
gets logT= 6.56i£^ (Table 3). 

The IUE satellite observed PKS 0558-504 twice with the 'Short Wavelength Primary' 
(SWP), low-resolution (AA = 6 A) spectrometer, on 1987 September 22, and 1989, November 
14-15. We retrieved the data products of these three observations from the MAST archive 1^1 
, and analyzed them with Sherpa (v. 4.1.2 ; Freeman, Doe & Siemiginowska, 2001). The 
three SWP Low-Resolution spectra (LR) of PKS 0558-504, are dominated by the intrinsic 
nuclear AGN continuum and emission line spectrum. However, the two spectra with highest 
S/N (obsid 37589 and 37604; Table 1) show also marginal evidence for a narrow absorption 
feature imprinted on the blue wing of the HI Lya broad emission line of PKS 0558-504, at 
A ~ 1360 A (Figure 10). We fitted the three spectra simultaneously with a model, consisting 
of a powerlaw plus narrow and broad components of the HI Lya and NV AGN emission lines 
and a negative absorption line to model the absorption feature at A ~ 1360 A. The best-fitting 
parameters for this absorption line are listed in Table 2. We tentatively identify this low- 
significance (1.7<t) absorber as an intervening broad HI Lya (BLA) at zbla = 0.119 ±0.001, 
a redshift consistent with that of the zblb HI absorber detected by FUSE. Figure 10 shows 
a portion of the highest S/N 1989 IUE-SWP LR (obsid 37604) spectrum of PKS 0558- 
504, with our best-fitting superimposed, and the residuals after subtracting the Gaussian at 
A = 1360 ± 1 A. 



4.1. Summary of the Observational Evidence 

The Far-UV and X-ray spectra of PKS 0558-504 show evidence for complex absoprtion, 
which is unlikely to be intrinsic to the AGN itself (§5.1). Figure 11 summarizes these 
conclusions, by showing the residuals (in a), in velocity space, to the best fitting continuum 
models to the FUSE (top panel), IUE-SWP LR (1989; second panel from the top) and 
XMM-Newton (GO: 3rd and 5th panels from the top; CAL: 4th panel from the top) spectra 
of PKS 0558-504. The residuals are centered around the absorption lines that we identify as 
HI hyp, HI Lya, OVIII Lya and FeXVII L (A = 15.015 A, rest frame), at a mean common 



12 http:/ /archive. stsci.edu/index. html 
13 http://cxc. harvard . edu/ciao/index.html 
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redshift of z — 0.118 ± 0.001 (all lines are consistent with this redshift). The red arrows in 
the figures indicate the best-fitting centroid of each line in their respective spectra. 

The FeXVII L A = 15.015 A) line (bottom panel) is only seen at 1.5a and 2a in the 
RGS1 GO and RGS2 CAL spectra (Figs. 5 and 6), and is not visible in either the RGS2 
GO or in the RGS1 CAL spectra, although all data are consistent with the presence of 
this line at the strength predicted by the best-fitting WHIM model. The OVIII Ka line is 
present both in the GO (3rd panel from the top) and CAL (4th panel from the top) RGS1 
spectra (the RGS2 is blind in this spectral region). Finally, possible HI lines are seen both in 
the FUSE (top panel) and IUE-SWP LR spectra: a BLB is detected in the high resolution 
FUSE spectra, at redshifts zblb = (0.1183 ± 0.0001), and an unresolved BLA is hinted in 
the IUE spectra at zbla — (0.119 ± 0.001). Both redshifts are consistent with that of the 
X-ray absorbers, within their la uncertainties. Finally we note that the measured ratio 
EWbla/EWbl£= (15 ± 10), between the equivalent width of the low-significance BLA line 
(IUE) and the BLB line (FUSE) is consistent, within the large statistical errors, with the 
predicted value of Hi)Lya/ hi)l v $— 7.4, for unsaturated lines . 

The combined statistical significance of this detection is 5.2a (4.6a if systematics are 
considered, and the small conribution from the putative HI BLA is not included). 

5. Discussion 

5.1. Ruling Out Intrinsic Absorption 

The Far-UV and X-ray absorber along the line of sight to PKS 0558-504, is only ~ —6000 
km s _1 from the systemic redshift of PKS 0558-504. In principle, then, it could be identified 
with a high-velocity outflow ejected from the Seyfert's nuclear region in the direction of 
our line of sight. Such phenomena are common in AGNs, and are known with the name 
of "Warm Absorbers (WAs, in X-rays: e.g. George et al. 1998; Piconcelli et al., 2005) 
or "Narrow Absorption Line" absorbers (NALs, in the UV: e.g. Crenshaw, Kraemer & 
George, 2003). However, we note that, the physical, chemical and dynamical properties of 
the detected absorber are all quite different from those typically found in WAs and NALs, 
commonly detected in Seyferts and higher luminosity quasars. 



1 For unsaturated lines, the equivalent width ratio of two electronic transitions 1 and 2 from the same 
series of the same ion X 1 , is linearly related to the product between their oscillator strength (/) ratio times 
the square of their rest-frame wavelength (A) ratio: i.e. (EWx*)i./(EWx«)2 — {fx i )i/(fx i )2 x (Ai/A 2 ) 2 . 
Thus, for the Lya and Lybeta transitions of HI, one expects: (EW Hi)hyal '(EW hi)l v i3— 7.4. 
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WA/NAL outflow velocities typically span a range of few hundreds to few thousand km 
s— 1, with >80% of the objects having t> OMt <1500 km s" 1 (e.g. Blustin et al., 2005 for WAs 
in X-rays, and Kriss, 2002 Crenshaw, Kraemer & George, 2003 for NALs in the UV) and are 
normally found to be present in several physical (temperature, and densities) and velocity 
components in a single object both in X-rays (e.g. Krongold et al., 2003, 2005, 2009) and in 
the UV (e.g. Kriss, 2002; Crenshaw, Kraemer & George, 2003). Here we measure an extreme 
velocity outflow of —6000 km s _1 and do not detect any other X-ray or FUV systems at lower 
velocity. The existence of very high ionization (FeXXV and FeXXVI) X-ray outflows, with 
sub-relativistic (i.e. ~ 0.1c) outflow velocities has been recently proposed for few objects 
(Reeves et al., 2009, and reference therein), but the existence and identification of these 
features are still highly debated (e.g. Uttley, 2009). In all such proposed cases, however, 
the high-ionization absorber has very high equivalent hydrogen column density (few x 10 23 
cm -2 ), is seen often together with other lower velocity, lower ionization components, and in 
some cases requires super-solar Fe abundance (Reeves et al., 2009 and references therein). We 
checked the high S/N XMM-Newton GO EPIC spectrum of PKS 0558-504 for the presence 
of blue-shifted He- or H-like Fe absorption, and did not find any down to an equivalent 
hydrogen column densities of >10 21 cm -2 . 

WA metallicities, when measured, are typically found to be at least Solar (e.g. Fields 
et al., 2005, 2007), in sharp contrast with the 1-5% Solar value we measure for this ab- 
sorber. We stress that our metallicity estimate for the z = 0.1 18 ±0.001 absorber is virtually 
model-independent. The HI fractions in purely photoionized or collisionally ionized clouds 
of gas that best-fit the RGS spectra of PKS 0558-504, are almost indistinguishable from 
the HI fraction that we derive from our best-fitting hybrid-ionization WHIM models. So, 
the metallicity-correction needed to reconcile the HI BLB (FUSE) and BLA (IUE) column 
density measurements with the X-ray equivalent hydrogen column density estimate, is inde- 
pendent on the model used to fit the X-ray data. 

Finally we stress that high-ionization WAs are virtually always found to co-exist with 
lower-ionization phases (e.g. HIP and LIP, Krongold et al., 2003, Andrade- Velasquez et al., 
2009), which in turn are often one to one associated with NAL components in the UV. In 
the case of PKS 0558-504 we do not detect any intrinsic OVI or CIV absorption (see also 
Dunn et al., 2007). All this makes the identification of the detected absorber with a high 
velocity nuclear ionized outflow, unlikely. 

A WA with unprecedented velocity, metallicity and ionization phases, appears to us to 
be less likely than the well-predicted WHIM possibility. We therefore tentatively identify 
this complex absorber with an intervening WHIM filament at the mean common redshift 
of z — 0.118 ± 0.001, and in the following discussion associate the X-ray absorber with the 
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Far-UV BLB and BLA detected in the FUSE and IUE spectra. 

5.2. Metallicity and Thickness of the WHIM Filament 

In principle, the detection of clearly associated metal and HI transitions from the same 
WHIM filament, allows one to measure the the metallicity of the absorber. In our case, 
the spectral resolution of the X-ray data is much poorer than that of the superior quality 
FUSE data, and the best-fitting width of the HI BLB absorber tentatively identified at a 
redshift consistent with that of the X-ray absorber (FWHM=160t|g km s -1 ) , is narrower 
than the expected HI BLB thermal width in gas with the temperature inferred from teh 
X-ray data. Thus a direct association between the HI BLB and the OVIII absorbers cannot 
be clearly established, based on these data. However, both the redshifts of the HI BLB and 
the OVIII absorbers, and the thermal width of the HI BLB lines, are consistent, within the 
large uncertainties (e.g. Fig. 9), with the hypothesis that the bulk of the HI and the OVIII 
absorbers are imprinted by the same gas. Under this assumption, we then estimate the 
metallicity of the putative WHIM filament, and note that this estimate would translate in a 
lower limit if the BLB and the metal absorbers were structured or not one-to-one associated. 

From the measured temperature, logT= 6. 56lo!o4 we derive the relative fraction of HI 
compared to total H: £ HI = (5.4l°;o) x As already mentioned in /S 3.3, the best- 

fitting temperature of the absorber is virtually independent on metallicity. We checked this 
by building two additional grids of models (other the one with Z = Z & used to fit the 
X-ray data) with metallicity of Z = 0.1Z Q and Z = O.O1Z . Figure 12 shows the OVIII 
to HI fractional abundance ratio as a function of the temperature of the gas, for Z = Z Q 
(black, solid curve), Z = 0.1Z Q (red, solid curve)and Z = O.O1Z (green, dashed curve). 
The ionization balance of the gas, at equilibrium, is virtually independent on the value of 
Z. From the metal-independent ionization correction value (i.e. £hi = (5-^-2x>) x 10~ 8 )> 
and the intrinsic (i.e. at z= 0.1183) HI column of N HI = (9 ± 2) x 10 13 cm -2 (infered 
from the measured EW(BLB)= 66 ± 16 mA and FWHM(BLB)= 1601^ km s" 1 ), we can 
thus derive the expected equivalent H column density of the absorber, which turns out 
to be: N H = (l.T^l) x 10 21 cm" 2 (Table 3). Our best fitting WHIM model measures 
N H = (3.2l?;£) x lO^Z/Z©)- 1 cm' 2 . This gives Z/Z Q = 0.02±°;° 2 , or Z = (1 - 4)% Z Q 
(Table 3). 

With an estimate of the equivalent hydrogen column density of the absorber in hand, we 
can estimate the thickness of the absorber along the line of sight. This is given, assuming ho- 
mogeneity by the ratio between the column and the volume baryon densities of the absorber: 
D =N b /n 6 . The value of n b is not constrained by the XMM-Newton data of PKS 0558-504. 
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However, hydrodynamical simulations that include feedback from galaxy super-winds predict 
that n b hi the WHIM correlates with temperature (e.g. Figure 3 in Cen & Ostriker, 2006), 
and that for a given temperature, lower metalicities are normally found in higher density 
environments (e.g. Figure 14 in Cen & Ostriker, 2006) . From Fig. 14 of Cen & Ostriker 
(2006), at T ~ 5x 10 6 K, metalicites of the order of \og{Z / Z Q ) ~ — 1.5 are reached at overden- 
sities, compared to the average density of the Universe (< n b >= 2x 10~ 7 (l + z) 3 (Q b h 2 /0.02)) 
of 5 = n b / < n b >~ 300 . At z = 0.118 and assuming h = 0.7, and f2 b = 0.046, this gives 
n b ~ 10~ 4 cm -3 . This implies D ~ 5^ x [n b / < n b >) _1 Mpc. 

Table 3 summarize the physical and geometrical parameters of the WHIM filament at 
z = 0.118 ± 0.001, along the line of sight to PKS 0558-504. 



5.3. Number Density and Cosmological Mass Density of WHIM 

A tentative OVIII detection along a single line of sight cannot be used to measure the 
number density of OVIII filaments and the cosmological mass density in the WHIM. However, 
we think that, in this particular case, this exercise is instructive, because it clearly shows the 
serendipitous nature of the detection, and its bias toward dense (i.e. high equivalent width) 
and rare absorbers. 

We estimate the number density of OVIII WHIM filaments, dN vin/dz and fi^ H/M , 
by taking into account the large Poissonian errors associated with small number statistics 
(Gehrels, 1986) . We get: dN OV iii/dz = 7.3^ 1 6 1 and Qf HIM = 0.18^ (90% sig- 
nificance). These numbers should be compared to the expected number density of WHIM 
filaments along a random line of sight, and to the total number of baryons as inferred by both 
microwave background anysotropies (Q b = 0.046 ± 0.002, Bennett et al., 2003; Spergel et al., 
2003) and by the standard "big-bang nucleosynthesis" when combined with light-element 



15 This is because super-winds are very efficient in metal-poiluting the IGM in the proximity of the host 
galaxy, but can only spread a finite amount of metals in a finite volume. Therefore IGM regions with lower 
primordial baryon densities will increase their gas metal content more than regions with initially higher 
baryon densities. 

16 This is consistent with our a-priory assumption of freezing the value of the baryon density to = 10~ 5 
cm~ 3 in our fit to the X-ray data, since at n& > 10~ 5 cm~ 3 the contribution of photoionization is negligible, 
and the best-fitting temperature and column densities are not affected by the exact value of this paremeter. 

17 We do not try to include systematics here, because of the aready big (and probably domaniting) statistical 
uncertainties associated with a single detection, which allow for the estimate to be consistent with all possible 
theoretical predictions. 
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ratios (Q b = 0.044 ± 0.004, Kirkman et al., 2003). From Cen & Fang (2006; their Figure 
4), we derive dNovin/dz(EW > 9 mA) ~ 0.8. This is consistent with the 90% significance 
lower limit of our estimate. Similarly, our estimate of Q^ HIM /Qf, = 3.9I35 greatly exceeds 
the expected value of about 40-50 % at z ~ (e.g. CO06), but is consistent within the 
large 90% uncertainities. This is not surprising, since the WHIM detection we report is a 
serendipitous discovery in the spectrum of a target whose observations were not planned with 
this aim, and as such is biased toward dense (i.e. high equivalent width) and rare absorbers. 

5.4. Caveats and their Solutions 

Although the association between UV and X-ray absorbers reported here, makes this 
detection of a probable WHIM filament the first for which a metallicity measure is possible, 
we caution that, due to instrumental limitations and uncertainties, a definitive identification 
of this absorber with an intervening WHIM filament, is not possible with the current data. 

First, the absorption that we identify here in X-ray as an OVIII Lya at z = 0.117±0.001, 
could also, in principle be identified as OVII K/3 near the systemic redshift of PKS 0558-504. 
We reject this interpretation based on the non-detection of associated OVII Ka absorption 
with the predicted strength ( > 6 times stronger than the detected line) - together with several 
other strong lines that should be detectable in the X-ray spectrum (CVI, NVI, NVII, NelX, 
etc.) - and the lack of intrinsic OVI absorption in the RGS and FUSE spectra of PKS 0558- 
504, respectively. Moreover, the absorption seen both in FUSE and IUE and here identified 
as HI absorption cannot plausibly be identified with any transition at or near the systemic 
redshift of the target. Unfortunately, however, the OVII Ka transition which would be 
associated with OVII K/3 at z ~ z sou , falls close to a strong RGS instrumental feature, which 
makes the source-frame absorption hypothesis difficult to reject with high confidence. The 
identification of the reported X-ray absorber with intrinsic nuclear OVII K/3 is thus unlikely, 
but not impossible based on the current X-ray data only. The only way of confirming 
the intervening origin of the X-ray absorber would be to clearly detect the other two lines 
predicted by the best-fitting WHIM model (one of which already hinted in the GO RGS1 
and CAL RGS2 spectra of PKS 0558-504): i.e. FeXVII L(A = 15.015 A) and NelX Ka, 
which fall at A = 15.02 and A = 16.77 A. This is a task that only the Chandra LETG grating 
could accomplish, thanks to its 40 cm 2 effective area at this energies (vs ~ 50 cm 2 for the 
RGS2 and ~ 15 cm 2 for the Chandra MEG), and factor of ~ 2.5 better spectral resolution, 
compared to the RGS. 

Second, while our identification of the absorber as a system at z — 0.118 ± 0.001 is 
the most probable with the currently available data, its origin might not be intervening. 
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The absorber could also be intrinsic to the nucelar Seyfert environment, outflowing from 
the nucleus of PKS 0558-504. We discuss this possibility in §5.1, and conclude that the 
physical parameters inferred for this absorber, make this interpretation unlikely. A way to 
definitely exclude this possibility, is to look for variability of the X-ray and far-UV features 
reported here. Observations with the HST-COS of the far-UV spectrum of PKS 0558- 
504, would not only definitely confirm or rule out the existence of Lya-(3 absorption at 
z = 0.1183 ± 0.0001 (the redshift in FUSE), but would also test for opacity (i.e. ionization 
state and/or column density) variability of this absorber, which, if detected, would make 
the WHIM identification impossible. Analogously, new high S/N X-ray (Chandra) high- 
resolution spectra of PKS 0558-504 would allow one to check for variability of the X-ray 
absorber. 

A third observation that can help sheding light on the true identification of this ab- 
sorber, by either strongly supporting or weakening the support for its intervening nature, 
is a measure of the galaxy density in the region surrounding the putative WHIM filament 
at z = 0.118. WHIM filaments are supposed to correlate strongly with concentrations of 
galaxies and with Large Scale Structure (e.g. Stocke et al., 2006b). Unfortunately, the line 
of sight to PKS 0558-504 lies in the Southern sky, at a location which is currently poorly 
investigated: none of the major galaxy surveys or deep surveys (e.g. 2dF , 6dF , SDSS 
), covers this region of the sky. Wide-field, multi-band observations, centered on the line 
of sight to PKS 0558-504, and sensitive down to a fraction of L* at z = 0.118, would allow 
photometric redshifts of the galaxies in the field to be determined with sufficient accuracy to 
exclude foreground and background objects. Spectroscopic follow-ups would then allow one 
to precisely estimate the galaxy density at z = 0.118 in a surrounding area around the line 
of sight to PKS 0558-504, to investigate on the association of WHIM filaments and LLSs. 



6. Conclusions 

We reported on the first combined and possibly associated Far-UV and X-ray tentative 
detection of a WHIM filament at a mean common redshift of z = 0.118 ± 0.001, along the 
line of sight to the Seyfert PKS 0558-504. Our main findings are: 

• OVIII Lya absorption is tentatively identified at z = 0.117 ±0.001 in two independent 



18 http: / /www. mso.anu.edu.au/2dFGRS / 
19 http://www. aao.gov.au/local/www/6df/ 
20 http:/ /www. sdss.org/ 
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XMM-Newton RGS1 spectra of the Seyfert galaxy PKS 0558-504: a high S/N 480 ks 
net GO spectrum taken in 2008, and the sum of a number of much lower S/N spectra 
taken between 2000-2001 as part of an XMM-Newton calibration campaign, with a total 
net expsoure of 309 ks. The combined, single-line significance of these two detections 
is 2.8(7. 

• When fitted with a model including proper parameterizations of the continua together 
with our hybrid (collisional ionization plus photoionization) WHIM absorption model, 
the two ful-band RGS1 and RGS2 spectra are consistent with the presence of a WHIM 
filament at z = 0.117 ±0.001 with logT= 6.56l{|;$ and N H = (3.2±^) x 1O 19 (Z/Z )- 1 
cm -2 . Based on X-ray data only, the parameters of this model are poorly constrained 
and both represent only lower limits at a 3a statistical confidence level. 

• Archival FUSE data of PKS 0558-504 show the presence of a broad absorption com- 
plex at A = 1147.1 ±0.1, with single-line statistical significance of 4.1a" (3.7<r when 
systematics are included). This line can be identified as broad HI Ly/3 at zblb = 
(0.1183 ±0.0001), consistent with the redshift of the putative X-ray OVIII Lya, within 
the large X-ray uncertainties. The width of this BLB absorber is marginally consis- 
tent, at a 2.5a level, with the thermal width of HI ions in gas with logT= 6.561q;{7 
(the best-fit temperature derived from the X-ray data). 

• Archival IUE-SWP, low-resolution, spectra of PKS 0558-504 taken in 1987 and 1989, 
also hint to the presence of an unresolved intervening HI Lja at z BLA = 0.119 ±0.001, 
with a single-line statistical significance of only 1.7a. The redshift of this absorber 
is consistent, within their relative la uncertainties, with the redshift zblb of the HI 
absorber detected with FUSE. The BLA (IUE) to BLB (FUSE) equivalent width ratio 
of these two lines, is consistent with the expected unsaturated value. 

• We tentatively associate the OVIII Lja X-ray absorber at z — 0.117 ± 0.001, with 
the far-UV BLB and BLA absorbers at z BLB = (0.1183 ± 0.0001) and z BLA = 0.119 ± 
0.001, respectively, and identify them with an intervening WHIM filament at the mean 
common redshift of z — 0.118 ± 0.001. The combined significance of this detection is 
5.2a (4.6a if systematics in the FUSE continuum modeling are taken into account and 
the IUE HI BLA line is not considered). 

• The above identification allows us to estimate for the first time the equivalent hydrogen 
column density of the absorber and so its metallicity under the assumption that the 
bulk of the HI BLB and the OVIII absorbers are physically associated. These turn out 
to be: Ntf = (1.51°;^) x 10 21 cm" 2 and Z = (1 - 4)%Z Q . If the HI BLB and the OVIII 
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absorbers are not directly associated or are structured, the above estimate represent 
upper and lower limits, respectively. 

• The non-detection of associated OVI absorption in the FUSE spectrum of PKS 0558- 
504, allows us to put a stringent lower limit on the temperature of the absorber: by 
combining X-ray and FUSE data we obtain: logT= (6.56+q;o4) K. 

• From the theoretical correlation between the temperature of WHIM filaments and 
their overdensity (relative to the average density in the Universe) and the expected 3D 
metallicity-temperature-overdensity relationship, we derive an overdensity of 5 ~ 300 
for our system, and therefore a thickness along the line of sight of 51^ Mpc. 

• Finally from this single detection we extrapolate the number density of OVIII Lya 
WHIM absorbers, and the cosmo logical mass density of WHIM. Both are consistent, 
within their large la uncertainties due to the low- number statistics, with the expected 
values, but their central-values are significantly higher due to the natural bias toward 
dense (i.e. high equivalent width) and rare absorbers associated with serendipitous 
WHIM detections. 
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Table 1: Log of the XMM.-Newton, FUSE and IUE Observations 



Dataset 


Date of Obs. 


Obs. ID 


Exposure" 


XMM-RGS CAL 


2000/02/07 


0116700301 fc 


19.0 




2000/02/10 


0117500201 


40.8 




2000/02/12 


0117710701 


49.9 




2000/02/13 


0117710601 


56.9 




2000/02/14 


0117710501 


6.8 




2000/03/07 


0120300501 


6.5 




2000/03/07 


0120300601 


7.0 




2000/03/07 


0120300801 


39.7 




2000/05/24 


0125110101 


30.0 




2000/10/10 


0129360201 


26.4 




2001/06/26 


0137550201 


14.8 




2001/10/19 


0137550601 


14.6 


XMM-RGS GO 


2008/09/07 


0555170201 


121.4 




2008/09/09 


0555170301 


127.8 




2008/09/11 


0555170401 


126.2 




2008/09/13 


0555170501 


126.6 




2008/09/13 


0555170601 


115.8 


FUSE-LWRS 


1999/12/10 


P1011504000 


46.2 




2002/11/07 


C1490601000 


48.3 


IUE-SWP LR 


1987/09/22 


SWP31899 


13.8 




1989/11/14 


SWP37589 


13.8 




1989/11/15 


SWP37604 


16.2 



a Net, in ksec. b RGS1 only. 
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Table 2: Absorption Lines (and upper limits) in the XMM-Newton RGS, FUSE 
and IUE spectra of PKS 0558-504 



Wavelength Width EW Significance 
in A FWHM in km s -1 in mA in a 


Id Rcdshift 

in km s _1 if z< 10~ 3 


XMM-Newton RGS 
23.52 ±0.03 < 900 23.7±|;? 6.4 
21.61 ±0.03 < 900 12.9 ±5.2 2.5 
21.17±0.03 < 1200 9.1 ±3.3 2.8 
15.00-15.03 a < 900 (fixed) < 3.0 l b 
16.76-16.79° < 900 (fixed) < 5.0 l b 


Ol Ka ±380 km s- 1 
OVll Ka 138 ±380 km s" 1 
OVIIIj{- Q 0.116 ±0.002 
NeIX Ka 0.117 ±0.001 
FeXVII L 0.117 ±0.001 


FUSE-LWRS 

1147.1 ±0.1 160±3Q 66 ±16 4.1 
1152.7-1154.7° 50 (fixed) < 12 l b 


BI Lyfj 0.1183 ±0.0001 
OVI(A = 1031.9) 0.118 ±0.001 


IUE-SWP LR 

1360 ±1 < 1300 1000 ±600 1.7 


}il Lya 0.119 ±0.001 


a Wavelength interval over which the la EW upper limit is computed, corresponding to z = 0.117 ± 0.001 
for the X-ray lines, and z = 0.118 ± 0.001 for the UV line. 



b la upper limit. 



Table 3: Physical and Geometrical Parameters of the WHIM Filament at z = 

0.118 ±0.001 



Rcdshift 


log(T) N ff Metallicity Thickness a 
(T in K) (in 10 21 cm" 2 ) in Mpc 


0.118 ±0.001 


6.56 +u 04 19 1.5l u i (1-5) %Z & 5l 2 



a For a baryon volume density of n& = 10 4 cm 
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20 21 22 23 24 

M}n A) 

Fig. 1. — 20-24 A portion of the RGS1 CAL (empty squares) and GO (filled circles) source 
(black and green , respectively) and background (red and blue, respectively) spectra of 
PKS 0558-504. The average background level during the CAL observations is about 10 
times higher than that of the GO observations. 
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X (in A) 



Fig. 2. — 20-24 A portion of the CAL (green points and errorbars) and GO (black points 
and errorbars) RGS1 spectra of PKS 0558-504, along with their best fitting continuum plus 
absorption models. The deep line-like features present in both spectra are instrumental 
features, due to bad-pixels/columns of the read out CCD. The remaining three real features 
present in both spectra are marked and here identified as 01 Ka and OVII Ka at z ~ 0, 
and, tentatively, as OVIII Lya at z = 0.116 ± 0.002. 
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14 16 18 20 22 24 

X(fn A) 

Fig. 3. — Best fitting WHIM model (black curve) together with its negative (red curve) and 
positive (blue curve) la temperature indeterminations. 
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Fig. 4. — 21-24 A portion of the RGS1 GO (top panel) and CAL (bottom panel) spectra of 
PKS 0558-504. The superimposed red curve is the best fitting WHIM model convolved with 
the instrumental responses. 
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Fig. 5. — 14-18 A portion of the RGS1 (top panel) and RGS2 (bottom panel) GO spectra 
of PKS 0558-504. The superimposed red curve is the best fitting WHIM model convolved 
with the instrumental responses. 
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Fig. 6. — 14-18 A portion of the RGS1 (top panel) and RGS2 (bottom panel) CAL spectra 
of PKS 0558-504. The superimposed red curve is the best fitting WHIM model convolved 
with the instrumental responses. 
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Fig. 7. — 68%, 90% and 95% logT - logN# contours for the putative X-ray WHIM filament 
at z = 0.117 ±0.001 
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Fig. 8.— 1137-1155 A portion of the FUSE spectrum of PKS 0558-504, in two different 
binning schemes: ~ 15 km s~l, the FUSE resolution (top panel), and ~ 90 km s _1 (bottom 
panel). Absorption lines identifications are labeled in both panels. Low- and high- velocity 
components of the low- ionization Fe and P and molecular H^absorbers are marked in blue 
and green and blue and cyan, respectively, but only the much stronger low- velocity lines are 
fitted (with the exception of the high-velocity FeII(A1144.94) line at A = 1145.7 A, which we 
fit because possibly partly contamined by intervening HI absorption). The broad absorption 
complex that we tentatively identify as intervening HI Ly/3 is clearly better seen in the 
heavily binned spectrum (bottom panel). 
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Fig. 9. — Gaussian normalization and FWHM 1, 2 and 3a contour plot, for the z 2 = 
(0.1183 ± 0.0001) HI Ly/3 absorber. FWHM=1.5 A corresponds to v = 390 km s" 1 at 
A = 1147.1 A, and so to the expected doppler parameter of b ~ 165 km s _1 , for hydrogen in 
gas at logT~ 6.5. 
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Fig. 10.— 1300-1500 A portion of the 1989, November 15 IUE-SWP LR spectrum of 
PKS 0558-504 (obsid 37604, Table 1) with our best-fitting superimposed (top panel), and 
the residuals after subtracting the Gaussian at A = 1360±1 A (bottom panel). Our tentative 
identification of a HI Lya absorber at z = 0.119 ± 0.001, is labeled. 
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Fig. 11. — Residuals (in a), in velocity space, to the best fitting continuum (plus ISM low- 
velocity lines, in the case of FUSE) models to the FUSE (first and second panel from the 
top), IUE-SW LR (1989; third panel from the top) and XMM-Newton (GO: 4rd and 6th 
panels from the top; CAL: 5th panel from the top) spectra of PKS 0558-504. The residuals 
are centered around the absorption lines that we identify as HI Ly/3, HI Lya, OVIII Lyo; 
and FeXVII L (A = 15.015 A, rest frame), at a mean common redshift of z = 0.118±0.001. 
The red arrows in the figures indicate the best-fitting centroid of each line in their respective 
spectra. The red segment labeled "OVL=o.ii8±o.ooi" in the second panel from the top, mark 
the spectral region where OVI associated with the HI BLB is expected, and not seen. The 
feature labeled "FPN" (Fixed Pattern Noise) in the first and second panel from the top, is 
a known FUSE Lif2A segment instrumental feature (e.g. Sembach et al., 2004) 
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Fig. 12. — OVIII to HI fractional abundance ratio (i.e. ionization correction) as a function 
of the temperature of the gas, for three different absolute metallicities: Z = Z & (black solid 
curve), Z = 0.1Z Q (red solid curve) and Z = O.O1Z (green dashed curve). The ionization 
balance of the gas, at equilibrium, is virtually independent on the value of Z, and so is the 
ionization correction. 



